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Abstract. In this paper, we construct a new family of random series defined on K , indexed 
by one scaling parameter and two Hurst-like exponents. The model is close to Takagi- 
Knopp functions, save for the fact that the underlying partitions of R D are not the usual 
dyadic meshes but random Voronoi tessellations generated by Poisson point processes. This 
approach leads us to a continuous function whose random graph is shown to be fractal with 
explicit and equal box and Hausdorff dimensions. The proof of this main result is based on 
several new distributional properties of the Poisson- Voronoi tessellation on the one hand, an 
estimate of the oscillations of the function coupled with an application of a Frostman-type 
lemma on the other hand. Finally, we introduce two related models and provide in particular 
a box-dimension calculation for a derived deterministic Takagi-Knopp series with hexagonal 
bases. 



Introduction 

The original Weierstrass series (see [38] ) is one of the most well-known examples of continuous 
but nowhere differentiable functions. Among the more general family of Weierstrass-type 
functions, the Takagi-Knopp series can be defined in one dimension as 

oo 

VxGM, K h {x) = ^22- nH A{2 n x) (1) 

n=0 

where A(x) = dist(x,Z) is the sawtooth -or pyramidal- function and H G (0, 1] is called the 
Hurst parameter of the signal. Introduced at the early beginning of the 20th century (see 
|34[ [21] ) , they have been extensively studied since then (see the two recent surveys [2j [23] ) . 

The construction of Kh is only based on two ingredients: a sequence of partitions of R (the 
dyadic meshes) associated with a decreasing sequence of amplitudes for the consecutive layers 
of pyramids. Therefore, we can easily extend definition ([1]) to dimension D ^ 2 using the 
Z?-dimensional dyadic meshes. 

In order to provide realistic models, it is needed to randomize such deterministic functions. 
Two common ways to do it are the following: either the pyramids are translated at each 
step by a random vector (see e.g. [3UJ ESI E])) or the height of each pyramid is randomly 
chosen (see in particular [15] for the famous construction of the Brownian bridge) . We obtain 
suitable models for highly irregular signals such as rough surfaces (see [151 123 E2] and [32] 
Chapter 6). 

In many cases, the graph of such functions is a fractal set. Therefore, their fractal dimensions 
constitute a crucial information for estimating the roughness of the data. The two most 
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common fractal dimensions are the box-dimension and the Hausdorff dimension. The former 
is in general easier to get whereas the latter is known only in very special cases (see e.g. 

[Ml [Ml EH 13 H). 

In this paper, a new family of Takagi-Knopp type series is introduced. Contrary to the previ- 
ous randomization procedures, our key-idea is to substitute a sequence of random partitions 
of M D for the dyadic meshes. An alternative idea would have been to keep the cubes and 
choose independently and uniformly in each cube each center of a pyramid. Notably because 
the mesh has only D directions, it would be very tricky to calculate a Hausdorff dimension. 
One advantage for applicational purposes may be to get rid of the rigid structure induced 
by the cubes and to provide more flexibility with the irregular pattern. A classical model of 
random partition is the Poisson-Voronoi tessellation. 

For a locally finite set of points called nuclei, we construct the associated Voronoi partition 
of MP by associating to each nucleus c its cell C c , i.e. the set of points which are closer to c 
than to any other nucleus. When the set of nuclei is a homogeneous Poisson point process, 
we speak of a Poisson-Voronoi tessellation [30} l2S~l [TJ. In particular, it is invariant under 
any transformation of W D which preserves the Lebesgue measure and its cells are almost 
surely convex polytopes. Classical results include limit theorems [I], distributional [51 [6] 
and asymptotic results [18^ [8] for the so-called typical Poisson-Voronoi cell. The model is 
commonly used in various domains such as molecular biology [31], thermal conductivity [22] 
or telecommunications (see e.g. [37] and Volume 1, Chapter 5 from [3]). The only parameter 
needed to describe the tessellation is the intensity A > 0, i.e. the mean number of nuclei 
or cells per unit volume. In particular, the mean area of a typical cell from the tessellation 
is A -1 . Multiplied by the scaling factor A"d, the Poisson-Voronoi tessellation of intensity 
A is equal in distribution to the Poisson-Voronoi tessellation of intensity one. This scaling 
invariance is a crucial property that will be widely used in the sequel. 

Let A > 1 and a, f3 > 0. The parameter A is roughly speaking a scaling factor and a, f3 
are Hurst-like exponents. For every n ^ 0, we denote by X n an homogeneous Poisson point 
process of intensity \ n/3 in IR D and by T n = {C c '■ c S X n } the set of cells of the underlying 

n/3 dcf r n /3 

Poisson-Voronoi tessellation. We recall that X~T n = \\~C C : c S X n \ is distributed as 7o 
and \~oT n = Tn-i thanks to the scaling invariance. 

Let A n : M D — > [0, 1] be the random pyramidal function satisfying A n = on UcgA 1 <^c 
and A n = 1 on X n (see Figure [Q . 

In particular we can then define a continuous function by 

oo 

V x E R D , F w (x) = A - ^ A n (x) . (2) 

n=0 

Let us denote by dime (A') and dimn(A) the (upper) box-dimension and the Hausdorff di- 
mension of a non-empty compact set K (see e.g. [TJ] for precise definitions). We are mainly 
interested in the exact values of these dimensions. Our result is the following: 

Theorem 1. Let A > 1 and < a ^ /3 ^ 1. Then F\ a g is a continuous function whose 
random graph 

Tx,a,p = {(x,Fx, a ,p(x)) :^[0,l] D }cK D xl 
is a fractal set satisfying almost surely 

ct 

dim B (I\ ajj g) = dim H (r A)a)( g) = D + 1 - -. (3) 
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(a) The Poisson-Voronoi tessellation T n . (b) The pyramidal function A 



Figure 1. Construction of the elementary piecewise linear function A n . 

Equalities © imply that the smaller % is, the more irregular F\^ a ^ and ^\ a ,p ar e (see Figure 
[2] and Figure [3]). The result of Theorem Q] naturally holds when [0, 1] D is replaced with any 
cube of R D . 
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(a) (A, a, /?) = (1.2, 1, 1) (b) (A, a, f3) = (1.2, 0.2, 1) 

Figure 2. Graph of the random function F\^ a ^ when D = 1. 

The paper is organized as follows. In the first section we state some preliminary results 
related to the geometry of the Poisson-Voronoi tessellations. We introduce in particular the 
oscillation sets O n ^ (see ([6])) that are used to derive explicit distributional properties on 
the increments of A n and precise estimates on the increments of F\ )Cl ^. Section 2 is then 
devoted to the proof of Theorem [H i.e. an upper bound for dinTr^r^Q.,^) comes from the 
estimation of the oscillations of -Fa,o,/3 whereas a lower bound for dmi^iT \ )0i ^) is obtained 
via a Frostman-type lemma. Finally, we introduce and study in the last section two related 
models: a deterministic series based on an hexagonal mesh and a random series based on a 
perturbation of the dyadic mesh. 

In the sequel we will drop the indices A, a and f3 so that F = F\ t0l ^ and T = T\ a g. 
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(a) (A, a, P) = (1.5, 1, 1) (b) (A, a, p) = (1.5, 0.2, 1) 

Figure 3. Graph of the random function Fx :tx ,p when D = 2. 



1. Preliminary results 

1.1. Notations. 

We consider the metric space M> D , D ^ 1, endowed with the Euclidean norm || • ||. The closed 
ball with center x € MP and radius r > is denoted by B r (x). We write Vol(.A) for the 
Lebesgue measure of a Borel set A C R D . In particular kd = Vol(.Bi(0)). The unit sphere 
of M D is denoted by S D ~ 1 and (Jd-i will be the unnormalized area measure of S^ -1 . The 
area surface of S D ~ l is then gjd-i = f.D-i('§' D_1 )- Finally, for all s 0, the s-dimensional 
Hausdorff measure is H s . 

For all x, y G [0, 1] D and all n ^ let 

Z n (x,y) = \-^(A n (x)-A n (y)) (4) 

so that F(x) - F(y) = Y^=o Z n(x , y) , and 

oo 

S n (x,y) = z m{x,y) (5) 

m=0 

so that F{x) — F(y) = Z n (x,y) + S n (x,y). Notice that Z n (x,y) and Z m (x,y) are two 
independent random variables for m ^ n. In particular Z n (x, y) and S n (x, y) are independent. 

Finally, we fix H > [3. For all n 0, we set r n = A~~d~. 

1.2. Random oscillations sets. 

Remember that the function A n is locally linear. Any maximal set on which A n is linear is 
the convex hull Conv({c} U /) of the union of a nucleus c from X n and a hyperface / of the 
cell associated with c. The set S n of such simplices tessellates M. D . For all n, N ^ we define 
the random sets 

On,N = {x G [0, 1] D : all points of B TN (x) are on the same simplex of S n as x} (6) 
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and 



w N = pi o n , n . 



(7) 



n=N 

The first result states that these sets are not too 'small'. 
Proposition 1.1. 

(i) There exists a constant C > such that, for all x € [0, 1] D and all n,N^0, 

nB-NH 

nxto n , N )^c\-^. 

(ii) We have limjv->oc P(Vol(Wjv) > 0) = 1. 
Proof. 

(i) By invariance by translation of X n and T n , we notice that for every x £ M. D , 

F(0?O ntN )=F(x<?O ntN ). (8) 

Let Sk n be the skeleton of the simplex tessellation S n , i.e. the union of the boundaries of all 
simplices (see the grey region on Figured]). 




Figure 4. The skeleton of the complete tessellation (in grey). 

In particular, we have the equivalence 

x^O n)N ^xeSk n +B TN (0) . (9) 

Let U be a uniform point in [0,1]'°, independent of the tessellation T n , and P# n be the 
distribution of the Poisson point process X n . Using ([8]) and Fubini's theorem, we get 

F(U£O n , N ) = [ F Xn (x<£O n>N )dx = F(0<£O n , N ). 
J[o,i] D 

Moreover, the equivalence ([9]) implies that 
P(0 O n , N ) = E Xn (^J ^ l Skn +BTN{0) (x)dx^j = E Xn (Vol((Sk n +B TN (0)) n [0, 1] D )) . (10) 
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It remains to calculate the Lebesgue measure of the set (Sk„ +B TN (0)) n [0, l] D . Denoting by 
J- n the set of hyperfaces of the simplex tessellation S n which intersect [0, 1] we have 

Vol((Sk n +B TN (0)) n [0, 1] D ) < Yo1 ^ + B th(°)) n [°. + Vol(a([0, 1] D ) + B TN (0)) 



< 2A-^ (2D + £ n [0, . 



Consequently, we have 

E* n (Vol((Sk n+J B TJV (0)) n [0,1] D )) < 2A~^ (E Xn (H D -\Sk n n[0,l) D ))+2D) . (11) 
Using the invariance of X n by scaling transformations and translations, we observe that 

E^(^-^Sk n n[o,i]°))=A- : ^E^(^- 1 (Skon[o,A^] B ); 

= A- 11 ^ E (W^-HSko n[0, 1] D )) (X^) D . (12) 

Combining (|10p . (jlip and (|12p . we obtain the required result (i). 
(ii) The point (i) implies that 

E(Vol([0, 1] D \ O n , n )) = f ¥(x^ O n , n )dx = P(0 £ O n , n ) = 0(A 1 ^). 

J[o,i] D 

Therefore, we obtain 

oo oo 

E(Vol([0, 1] D \ W N )) ^ E ( Vol ([0, 1] D \ On,n)) < X) °( a! ^) = CKA^^). 

n=N n=N 

Finally, using Markov's inequality, 

P(Vol(Wjv) < 1/2) = P(Vol([0, 1} D \ W N ) > 1/2) < 2E(Vol([0, 1] D \ W N )) < 0(A^^) 
and the result (ii) follows. □ 



I. 3. Distribution of the random variable Z n (x,y). 

This subsection is devoted to the calculation of the distribution of Z n (x,y) conditionally 
on {x £ O n ,n} when \\x — y\\ ^ r n . In particular, we obtain in Proposition 11.21 below 
an explicit formula and an upper-bound for the conditional density of Z n (x,y). A similar 
method provides in Proposition II. 31 the integrability of the local Lipschitz constant of Zq. All 
these results will play a major role in the estimation of the oscillations of F (see Proposition 

II. 4p and the application of the Frostman criterion (see Proposition ll.5p . 

For any x 6 [0, 1) D , let c n (x) (resp. C n (x)) be the nucleus (resp. the cell) from the Voronoi 
tessellation T n associated with x. Let c' n (x) be the 'secondary nucleus' of x, i.e. the point 
of X n which is the nucleus of the neighboring cell of C n {x) in the direction of the half- line 
[cn(x),x). Moreover, for any z\ / Z2 € M D and x G M D \ {z±, Z2}, we consider 

- the bisecting hyperplane H Zl ,z 2 of [zi,Z2], 

- the cone A(z±,x) centered at z\ and generated by the ball B Tn (x), 
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the set A n>x of couples {z\ , z 2 ) with z\ G" B Tn (x) such that x is between the hyperplane 
orthogonal to z 2 — Z\ and containing z\ and the parallel hyperplane which is at distance 
r n from -Hzi,z 2 on the zi-side: 



A n , x = {( Zl ,z 2 ) £ B Tn (x) c xR D : < (x - z u z 2 - Zl )^ \ \\z 2 - zi|| 2 (l - .. 2? " n „ ) \. 

I 1 ll z 2 — J 

Finally, we denote by V n (x, Z\, z 2 ) the volume of the Voronoi flower associated with the 
intersection H zljZ2 n A(zx, x): 

V n (x,zi,z 2 ) = Vol ( (J : u G H» 1)2a n A(^i,x] 




Figure 5. The configuration of (21,22) with the associated Voronoi flower (in red). 



Proposition 1.2. Let n ^ and cc,y G [0, 1] D suc/i i/iai x G On,n and < ||x — y\\ ^ r n . 
T/ien 

(i) The increment Z n (x,y) is given by 

Z n (x,y) = 71 ~( \ r\u2 ( X ~ y> C 'n( X ) ~ c n{x)). (13) 

(ii) The density z n of Z n (x,y) conditionally on {x G O n ,n} is given by ()17|) /or D ^ 2 
and by f)16|) /or D = 1. Moreover, it satisfies 

snpz n (t)^ J x-y\\ L X ° (14) 

where C is a positive constant depending only on dimension D. 

Proof. 

(i) If x G O njn and ||x — y\\ ^ r n then 

dist(x,g Cn ( a .) ;C / i (. r )) 
dist{c n (x), H Cn{x) ^ x) ) 

Moreover, 

+ c4(x) c n (x) - d n (x) 



A n {x) 

dist(x,H Cn{x)A(x) ) = 



||Cn(x) Cn(^) 
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and 

dist^^),^^)^^,)) = - \\c n (x) - c' n (x)\\. 
It remains to use the definition Q of Z n (x,y) to obtain the result (i). 



e„(x) V 




Figure 6. The random variable Z n (x,y). 

(ii) We need to determine the joint distribution of (c n (x), d n (x)). We notice that x belongs 
to On t n if and only if all the points of B Tn (x) have the same nucleus as x in X n and same 
'secondary nucleus'. In other words, 

x G O ntn <=>- A(c n (x),x) D flc»(aO,<(a:) c C Cn ^ H C^ x y 

By Mecke-Slivnyak's formula (see Corollary 3.2.3 in |33j). for any measurable and non- 
negative function h : M? — > R + we have 



E(h(c n (x),d n (x))l {x€0nn} ) 

= E ( S ^( 2; l' Z 2)l(c Tl (x),<(x))(^1^2)l{ a;6 O njn } 



' 21^22 



X2n jJ h ( z i> z *) ¥ ( A -( z i> x ) nH *i>* ^C Zi nC Z2 )l An Jz 1 ,z 2 )dz 1 dz 2 
\ 2np If h(z 1 ,z 2 )e W (-X n ^V n (x,z 1 ,Z2))lA n J z ^ z ^zidz 2 . 



We point out a small abuse of notation above: the sets C Zl and C Z2 are Voronoi cells associated 
with z\ and z 2 when the underlying set of nuclei is X n U {zi,z 2 }. We proceed now with the 
change of variables z 2 = z\ + pu with p > and u G S 15-1 . 

Let L n (x) = \\c' n (x) - Cn(x)\\, u n (x) = c " ( ^~^ (x) and V^(-) = V n (x, zi, zi + •) (keep in mind 
that V' n will still depend on x and z\ though this dependency will not be visible for sake of 
readability). The density of (L n (x),u n (x)) conditionally on {x G O n ,n} with respect to the 
product measure dpda£)_i(u) is 

prT(=f? \ / e^{-X^VM)l l o 4 - Tn] ((x-z 1 ,u))p D - 1 dz 1 . 
r{x g u n)Tlj ) Jb t (x) c 
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Using CGI), we can rewrite the quantity Z n (x,y) in function of L n (x) and u n {x) as 

*<,,,>--£»<,-,,..<,>>. 

The density of the distribution of Z n (x, y) conditionally on {x G O n ,n} can then be calculated 
in the following way: for any non- negative measurable function ip : R+ — > M+, 

A 2n/3 , 

K(ip(Z n (x,y))) = — — -/ J n (x,y,zi,p,n)dzi (15) 

r(x g u„ in j J B Tn (x)° 

where J n (x, y, zi, p, u) is equal to 



^(-2A-^p" 1 (x - y,u))e- xn0v ^ u H R+ ((x - z^p^dpdan-iiu), 
the domain of integration for p being [2r n + 2(x — zi, u), oo). 

Case D = 1. We observe that u = ±1 and V^(pw) = P- If u = 1 (resp. u = — 1), the range of 
zi is (— oo, x — r n ) (resp. (x + r n , oo)) while for fixed z\ the range of p is (2r n + 2(x — zi), oo) 
(resp. (2r n + 2(zi — x),oo)). Consequently, we have 

/ ^{2X—p-\x - y))e-»*\* ( M - 2r n ) dp. 



E(^(Z n (x,y))) 



P(X G O n ,n) 7|p|>4r„ 

Applying the change of variables £ = 2X~ na p~ 1 (x — y), we get that the density of Z n (x, y) is, 
for \t\ < \~ nal -^, 

2X 2np 2X n V - a) ]^vj ( \x- y\ \ na \x-y\ 

Zn(t) - F(x i O n , n f T" 2Tn J A ■ (16) 

In particular, 

^nW < =7 ^ r ; — SUp(e _2n « 3 ), 

P(x G C n ,n) «>0 

which shows fl" 



Case D ^ 2. We go back to (fT5|) . For almost any u G S D 1 , there exists a unique 
v G n {y — x} 1 - and a unique s = (u, jT^Efir )g ( — 1> 1) such that 

y-x r, 2 

U = U SV = St. + V 1 — V. 

\\y - x \\ 

In particular, we can rewrite the uniform measure of S 15-1 as 

dcr D _i(u) = (1 - s 2 )^dsrfcr£i_2(u). 
We thus get that J n (x, y, z\, p, u) is also equal to 

V>(2A-^||x-y|| S p-V An " K(pUs '^^ 



the domain of integration for p being [2r n + 2(x — zi,u S)V ),oo). 

We now proceed with the change of variables p = 2A~~d~||x — with st > 0. We then 

deduce that the density z n (t) at point t of Z n (x,y) conditionally on {x G O n ,n} is given by 

Z ^ = P( X ^ O ) j j j ^n^'V' S ' V ) 1 D n (x,y, Z!,t, S,v)dsda D -2(v)dZl ( 17 ) 
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where 

J>(*, y, ,1 , t, s, v) = e -** * (^%^) ( 2A ^";-^) " 



and 



A -C x — y \ \s 

D n = <j (x,y, z±,t, s,v) : \\x - zi\\ > r n and ^ (x - zi,u s>v ) < t t 



In the sequel, we only deal with the case t > but the same could be done likewise for t < 0. 
We denote by x' the intersection of the half-line [z±, x) with the boundary of the Voronoi cell 
of Z\. Moreover, we write z\ = x — where 7 > r n and w £ S D_1 . In particular, we notice 
that 

11 / P A o ||x — y\\s 

2(w,u s , v ) t(w,u SjV ) 

We can now easily estimate the volume V' n (-) in the following way: 

_ I, \ D 

KM > voi,iV-„, - "*( (< JXH ■ <■*> 

We then proceed with the following change of variables: for almost any w £ S D_1 , there exist 
a unique £ = (w, u S)V ) £ [0, 1) and a unique 7/ £ S D ~ 1 n{n Sj t,}- L such that it! = £« s ,«+"\/l — £ 2 r] 
and 

dz! = 7 D ~ 1 d7 dcr = 7 D ~ 1 (1 - £ 2 )^d 7 df d<7 D _ 2 (»7). (19) 
In particular, when (x, y, z\, t, s, v) £ D n , we have 

(x — Zi.Usv) A~~d~ \\x — y\\s , , 

^ 7= \ i^w^ 11 y\\ _ (2Q) 

Consequently, for fixed s,£ £ (0, 1), we have 

2 / \ _S2L 11 || \ ^ 

lz )n d7da D _ 2 (^)da D _ 2 (r ? ) < ^ / j (21) 

We deduce from (fT7|) . (fTBj) . ([THD and ([2"T]) that the density z n (i) satisfies, for every t > 0, 

A 2n^ w 2 ,1 ,1 




z n (t)^ nWrra D ~ 2 , / 4{x,y,t,a,0d8dt (22) 



where 



2-D B_ 3 



We then use the change of variables r = A~~d~||:e — i/||s(t£) _1 . 



2 
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Subcase D ^ 3. Using that (1 — s 2 )(l — £ 2 ) and £ are bounded by 1, we obtain from the 
change of variables that 



Znit) < 



X 2n/3 



D-2 



DF(x g o n>n ) 



D 



J{t>0} 



t x-^f 



\x - y\ 



dr d£ 



C 



(x G O n>n ) ||x - 2/|| 7{ T>0 } 



■ n(/3-q) 
A D 



P(xGO n ,„) ||x-y|| 
where C and C" are two positive constants which only depend on D. 

Subcase D = 2. We return to ([22]) and apply the same change of variables r = A~~d~||x 
y\\s(t£) . We now obtain that 

Zn(t) ^ 



where 



(...) 



1A 



x - y\\F{x G O n>n ) J {T>0} 



r V A ^(...)dr 



(23) 



— n f* 
A ~2~ \\x — y 



5=0 



A 2 ||x — y\ 



We notice that there exists a positive constant C > such that for every a > 0, we have 



lAa / t> 

(l-£)-5 l-£) d£<C|log|a-l||. 



(24) 



Indeed, due to the facts that the left-hand side of ([24"]) is bounded for large a and that the 
calculation is symmetric with respect to 1, it is enough to look for the behaviour of the 
Abelian-type integral when a > 1 is close to 1. A direct calculation shows then that 



(1 -£)-' ( 1 - -) 2 d£ = v / aargch 



a + 1 



a — 1 J a-+l 

which proves (]24p . Consequently, we get from (j23|) and (|24p that 



~ — log (a — 1) 



Zn(t) < 



|x - y||P(x G O n , n ) J {T>0} 

where C denotes again a positive constant. 
We now fix e G (0, 1) and we split the integral: 



log 



1 



A 2 ||x — y\ 



tr 



dr, 



on the range of r which satisfy |1 \\ x -y\\ | -> £) upper-bound is similar to the case 



D > 3; 

- on the range of r satisfying 1 1 



1-e) 



- e (l + e) 2 * 2 



A 2 | x — y| 



i < e, the integral is bounded by 

5 



1+e 



U=l — E 



log 1 1 - u-^du CV 



A 2 |sc — y\ 



where <p(u) = e " (!+ e ) 2 u 5 , it > 0. 
It remains to notice that the maximum of the function ip is of order 0(A~2 n ^) to deduce the 
required result (TH]) . □ 
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We conclude this subsection with the integrability of the Lipschitz constant L(x) of the affine 
part of Ao above x, i.e. 

L{X) = l|co(x)-c' (x)||- (25) 
Proposition 1.3. For every x G [0, 1] D , E(L(x)) < oo and sup E(L(x)| x G 0n,7v) < oo. 

n,7VeN 

Proof. We start by noticing that the conditioning can be easily overlooked. Indeed, the 
condition {0 G O n ^} implies that a vicinity of the origin is in the same simplex of Sq so 
co(0) and c' (0) are farther. Consequently, we have E(L(x) \ x G O n ,Af) ^ E(L(x)) for every 
n, N ^ 0. 

For D = 1, the integrability of the variable given by (|25p comes from the fact that the 
distance from the two neighbors of x (the nearest and second nearest) is Gamma-distributed. 
When D ^ 2, we use a reasoning similar to the proof of Proposition 11.21 to obtain that 



E(L(z))=E^ 1 Z2 || JMs),^))^!' *z) J 



■ — - — -P((x + R+(x - n H ZuZ2 G C 2l n C Z2 )d Zl dz 2 . (26) 

PI — Z2\\ 

We write z 2 = z\+pu where u G Sr 0-1 , p > and u = s + "^1 ~ s2t> where s G (0, 1) and 

v G §' D ~ 1 n{a; — zi} -1 ". In particular, the distance from z\ to the point (x+M+(x — zi))nH ZljZ2 
is . Consequently, we deduce from a change of variables applied to the integral in (|26p that 

^j^-1r + ( P - 2\\x\\s)e- KD (£) D p D ~\i - s^^dsdpdx. 

When D 3, we proceed with the change of variables r = There is a constant C > 
such that 



E(L(x)) 




oo 



E(L(x)) < o; D _2 / p " 1 / (/ dx)e-^ M ^ dp 



p=0 7r=| \JB{0,t) 



2 

OO /'OO 



,D dT 



r- 



< C I p D - x I T D - 2 e- KDT ° drdp 
'p=0 



-E 
' 2 

(•OO 



<cf ^e-TFCf )"dp f° r^V^dr < oo. 

Jp=0 ir=0 



lp=0 

Finally, when D = 2, with the same change of variables, we get 



E(L) <c r,r.-*-^=. (27) 

^=1 A /i_(^) 2 



(£)' 

We treat separately the integral in r for fixed p > : 



_2 



~ 2 



arccos 



P_ 

2u 



, 3 



OG 



T = 



< 2pe~ 7r U 

^^ pe -(f) 2 +C7 'e-^ ; (28) 
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where C is a positive constant. Inserting ([28)) in (|2"7|). we get the required result. □ 
1.4. Size of the increments of F. 

The box-dimension of T, as well as its Hausdorff dimension, is closely related to the oscillations 
of F (see [EH E]). Let us recall that, for every A C [0, 1]"°, the oscillation of F over A is 
defined by 

osc(A) = sup \F(y') - F(y)\. (29) 

y,y'€A 

In particular, we will consider, for all x G [0, 1] 15 and r > 0, the oscillation of -F over the cube 
x + [0, t] d given by 

osc r (x) = osc(x+ [0,r] D ) = sup \F(y') - F(y)\. (30) 

ifjj/'eaH-pjTp 

Proposition 1.4. Let < p < jj. Then, for every x G [0, 1] 15 , we Ziaue, w/ien TV — >■ oo, 

E(osc rjv (x)) = E(osc Tjv (0)) = 0(rP ). 

Proof. Let us write 

= sup \F(x) -F{0)\. 

We claim that 

P(liminf {S N ^ t^}) = 1. (31) 

By Markov's inequality, 

V{5n>t p n ) ^T- p E(5 N ). (32) 

We can write 

E(5 N )^E(jr sup \Z n (x,0)\)=S 1 (N) + S 2 (N) + S 3 (N) 

\ n=0 x<E[0,T N ] D J 

with 



S 1 (iV) = ^E( / sup |Z n (s,0)| G O n ,ivV(0 G C»„,jv), 
S 2 (N) = J2^( sup |Z n (x,0)| O n ,ivV(0 O n ,jv), 
S 3 (N)= e( sup |Z n (s,0)|Y 

n=N+l W[0,tjv] d / 

Since [0, T n ] D C By/j) T (0), we mention that for the purpose of this proof, the definition ([6]) of 

the set On,N should be slightly adapted by substituting \^Dtn for t^. For sake of simplicity, 
we omit that technical detail. 
For Si(N), we notice that 

G O ntN ) {\-^VDT N ){\t E(L(0) I G n ,jv)) 



E sup |Z„(x,0)| 
Vselo.TAr] 15 

where L(0) is the Lipschitz constant of Aq at when the underlying Poisson point process 
is homogeneous of intensity 1. Thanks to Proposition II. 3[ we have sup n jyeN E(L(0) | G 
O n M) < oo. 
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Thus, using < a ^ j3 ^ 1, we obtain 

N ( (P-a-H)N 

, . -r-^ na NH n/3 , Cl A D if CK < 5 

5i(JV)< VA"tA~ At Vd sup E(L(0) | e O n iv) < < —nh 

~^ n,N&N [ C[ NX^~ if a = (3 

where C\, C[ are two positive constants which do not depend on N. 
For 5^2 (AT") we use the upper estimate (see (|1D) 



E( sup \Z n (x,i 



and Proposition ll.lf i) to get 

V — "V na np— iV-H 



na ^ nl3-NH | C2 A D if « < /? 

— iVJf 

C 2 N\-D~ ii a = /3 



n=0 



where C2, C 2 are two positive constants which do not depend on N. 
Finally, for Ss(N) we only use the upper estimate 



E sup |4(z,0)|KA~ 

\x£[0,tn] d J 

to get 

00 

n=JV+l 

where C3 > is a positive constant which does not depend on N. 
Therefore, 

E > sup |Z n (x,0)| U < \ N N> « 

xe[o,r N ]B J [ C'{\log(T N )\T N + T») ifa = p 

where C, C are two positive constants which do not depend on N. The rhs of (132[) is in 
particular summable in N as soon as p < jj (which guarantees that 1 — > p since 
1 — -Tj > 0). Consequently, by Borel-Cantelli's lemma, (|3ip holds. Then, we obtain 

osc TiV (0) ^25 n = 0{t p n ) 

almost surely for all TV" large enough. To conclude, let us notice that, by stationarity, it 
follows that 

E(osc Tjv (x)) = E(osc TjV (0)) = 0(t p n ). 
for all x G [0, 1] D . □ 

We conclude this subsection with an estimate of the expectation of a particular functional of 
the increment F(x) — F(y) that will appear in the application of a Frostman-type lemma in 
the next section. 

Proposition 1.5. Let s > 1 and n 0. If x,y G [0, 1] D satisfy r n +i < \\x — y\\ ^ r n i/ien 

E((|F(x)-F(y)| 2 + ||x-2/|| 2 )-fl 0nin ( 2 ;)) < C\\x - y\\- s+ ^r (33) 
where C > is a constant which do not depend on x, y nor on n. 
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Proof. Remember that F(x) — F{y) = Z n (x,y) + S n (x,y) where Z n (x,y) and S n (x,y) are 
independent. Let P n be the probability associated with X n and fj,s n be the probability 
distribution of the random variable S n (x,y). From Proposition 11.21 one obtains 

E ((\F(x) - F(y)\ 2 + \\x - y\\ 2 )-^o n Jx) 

((Z n (x,y) + v) 2 + ||z - y\\ 2 yi H Gn n (x)dP n d/i 5n («) 
(x E O n>n ){(u + v) 2 + \\x - y\\ 2 )-h n (u)dud/is n (v) 



^ fp(i€ O n , n ) f / 



||x-y|| s z n (ii)ditd/i5 n (t;) + / + s z n (u)dudfi Sn (v) 

{|u+u|<||a;-y||} J {|w+«|>||a!-y||} 

^ 2\\x - y\\ sup(z n (f))||x - y\\~ s + sup(z n (i)) / \u + v\~ a dudfj,s n (v) 

t£M ieR ./{|u+u|>||a:-2/||} 

^ n n 1 n(ff-o=) I - i w(/3-oQ ii — .J-I 

< 2C ||x — 2/ 1 1 \\x — y\\ A d ||x — y|| + G||x-y|| A o ||x — y|| ^ 

n(/3-g) 

= C||x - 7/ 1| ~ s A b . 

Finally, the assumption on ||x — y|| implies that A _n ^ (||x — y||A^)», which provides the 
desired bound. □ 



2. Proof of the main theorem 

Let us recall that for any non-empty compact set K one has ^ dimn(A') ^ diniB(-fQ (see 
|14j). Thus the proof of Theorem [T] will consist in proving that D + 1 — % is an upper bound 
for dimB(r) and a lower bound for dimn(r). 



2.1. An upper bound for the box-dimension of T. 

First we investigate the box-dimension of T. For every r > we cover [0, 1] D x R with r-mesh 
cubes and denote by M(t) the (finite) number of cubes from this partition which intersect T. 
Then, we can express the box-dimension of T in terms of jV(tjv) (see Section 3.1 in [T3] and 
Section 2.2 in |36j): 

a- tr\ v logAT(rjv) ,„.s 
dime (I J = hmsup — - . (34) 

AT^oo -log (tat) 



Lemma 2.1. There exists a constant d > such that P(dimB(r) = d) = 1. 

Proof. For every m E N, let us denote by A m the cx-algebra generated by the point process 
X m . We observe that the graph of the function which associates to any x E [0, 1] D the quantity 
F(x) — Yl™=o ^ w &-n{x) has the same box-dimension as the graph of F. Consequently, we 
can use (|34p to show that dimB(r) is a random variable which is measurable with respect to 
m : m n) for every m E N. We then use the 0-1 law to deduce that it is almost surely 
constant. □ 

Proof of the upper bound of ([3|). 

Using (f34l) as if limsup^r^oo is a real limit, provided that tn is replaced with a subsequence, 
applying first Lebesgue's convergence theorem to ^ D + 1 and finally Jensen's 
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inequality to the concave function log, we obtain 

E(dim B (I>))= lim e(^P4)^ lim log if UM 



7V->oo \— log(TAr) / 7V->oo \— log(rjv). 

We now use Proposition 11.1 from [14] and Proposition 11.41 to get 

E(Af(r N )) < 2\t^] d + t n 1 Yl E(osc rjv (krAr)) = O^" 1 ). 

k=(fci,. ..,k D ) 

In conclusion, for every H > j3 and p < %, we have E(diniB(r)) ^ D + 1 — p. We finally get 

dim B (r) = E(dim B (r)) < D + 1 - ^ 
by letting H — > /3, p — > %, and by using Lemma l2.1i □ 

2.2. A lower bound for the Hausdorff dimension of T. 

To find a lower bound for the Hausdorff dimension of a compact set is generally a difficult 
problem. An important step was made in [19] when Hunt proposed a way to find a lower 
bound for the Hausdorff dimension of the graph of Weierstrass functions using the 'finite 
energy criterion'. This well-known criterion is used for calculating the Hausdorff dimension 
of more general fractal sets (see e.g. [29], Chapter 4). However, the arguments of [19] cannot 
be applied for all Weierstrass- type functions, in particular it rules out the Takagi-Knopp 
series because the sawtooth function A is not regular enough. 

Let us recall that the Hausdorff dimension of a non-empty compact set K may be expressed 
in terms of finite energy of some measures thanks to a lemma due to Frostman (see e.g. 
[201 [25]): 

dimn(-ft') = sup (sup{s ^ : J s (/i) < oo}) 

where the supremum on \i is taken over all the finite and non-null Borel measures such that 
fi(K) > 0, is(ju) being the s-energy of p, defined by 



I s {p) =JJ\\x- y\\- s d»(x)dfi(y). (35) 
Therefore, if such a measure \i satisfies < oo then dimn(i^) ^ s. 

Let us recall now a classical way to construct such a measure on V. For each N ^ the 
set Wn is a Borel subset of [0, 1]"°. Since F is a continuous function then T is a Borel set 
too so that we can consider the measure \iw N obtained by lifting onto T the D-dimensional 
Lebesgue measure restricted to Wn- Precisely, for all Borel set E C 

H Wn (E) = Vol{x G [0, 1] D n W N such that (x,F(x)) G E) 

and nw N is a positive measure as soon as Vol(Wjv) > 0. Using Euclidean norm on M. D x R 
we can recast (|35|) into 

I s {vw N )=ff {\F(x)-F(y)\ 2 + \\x-y\\ 2 )-Mxdy. (36) 

J J W N x W N 

Notice that the finiteness of I s (^w N ) depends only on the size of the increments F(x) — F(y) 
when \\x — y\\ is small. 
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Proof of the lower bound of ([3]). 

Let us consider, for all n ^ 0, the set 

T n = {(x,y) e [0,1] D x [0,1] D : r n+ i < ^ r n } . 

Let N ^ 1 and s > 1. Since Wn C ra ,n for all n ^ N, we have, with (j36| 

J-(^) < // i\F(x) - F(y)\ 2 + \\x - y\\ 2 )-^ dxdy 

J J(x,y)&W N x[0,l] D 



n= ■/ ■/ (x l3/ ) 6T„n ( Wjv x [0, 1] D ) 



// {\F(x)-F(y)\ 2 + \\x-y\\ 2 ) > l 0n Jx) dxdy 

where 

C N = ff (\F(x) - F(y)\ 2 + ||a; - y\\ 2 Y^ dxdy 

J J(x, ?/ )G(ToU-UT JV _ 1 )n(W/ J vx[0,l] c ) 

< // Ik - yH~ s dxdy < X NHs . 



JJix^eiToU-UTN-!) 

To show that the integral (|36p is finite almost surely it is enough to show that its expectation 
is finite. By Fubini's theorem, 

m^w N )) < \ NH ' + E // E ((l^) - F (y)l 2 + II* - y|lT f W*))^- 

n=N J J(x,y)eT n 
Using the maj oration (j33j) we obtain 



E(J fl ( m J) < \"» s + C \\x- y\r s+ — dxdy. 

J J (x,y)eu™ =N T n 

This latter integral converges as soon as — s + > —D. Therefore the random measure 
Hw N has a finite s-energy for all 1 < s < D + Since N may be chosen such that 

Vol (Wat) > with an arbitrary high probability (see Proposition II. ip we deduce that 

dim H (r)^ J D + ^^ 

for all H > P and almost surely. We obtain the desired lower bound by letting H £ Q + go 
to p. □ 



3. Related models 

In this section, we study the fractal properties of two different models which are related to our 
Poisson-Voronoi construction: a deterministic series of pyramidal functions with hexagonal 
bases on the one hand, a random perturbation of the classical Takagi-Knopp series on a 
dyadic mesh on the other hand. 
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3.1. Takagi-like series directed by hexagonal Voronoi tessellations. 

The series that we study here is very close to the original function F\ a g. The novelty lies 
in the A n signals: we consider now pyramids with a regular hexagonal basis. This model is 
naturally related to the previous one for two reasons: first an hexagonal mesh is classically 
known to be the Voronoi tessellation generated by a regular triangular mesh. Secondly, the 
mean of the number of vertices of a typical cell from a Poisson- Voronoi tessellation is known 
to be 6 (see e.g. [28], Prop. 3.3.1.) so that an hexagonal mesh may be seen as an idealized 
realization of a Poisson- Voronoi tessellation. 

We start with the deterministic Voronoi tessellation whose cells are identical regular hexagons 
such that one is centered at the origin (0, 0) and has a vertex at (1,0). Then, considering all 
the centers of these hexagons as a point process Xq, we set X n = 1~ u Xq and construct the 
hexagonal Voronoi tessellation associated with. Here again A n : M? — > [0, 1] is the pyramidal 
function satisfying A n = on (Jce-*^ ^ c an< ^ ^ n = 1 on <^n- 

Let us notice that, for all i£l 2 and all n ^ 0, we have A n (x) = A(2 n x). We fix a G (0, 1] 
and define a function 

oo oo 

f a (x) = 2~ na A n (x) = 2- nQ A(2"x). (37) 

n=0 n=0 

The main theorem of this section is the analogue of Theorem [TJ Actually we state a more 
precise result than Proposition 11.41 for the oscillation of f a but we cannot determine the 
Hausdorff dimension of T a . 

Theorem 3.1. Let < a ^ 1. Then f a is a continuous function such that 

3C,C">0, Vt€(0,1), VxG[0,l] 2 , C't q < osc r (x) < Cr a . (38) 
Moreover, its random graph 

r Q = {(*, f a (x)) : x g [0, l] 2 } Cl 2 xl 

is a fractal set satisfying 

dim B (r Q ) = 3 - a. (39) 

Proof. In the sequel we drop again the index a so that f = f a and r = T a . We also keep the 
notation Z n (x,y) = 2" nQ (A n (x) - A n (y)) = 2" nQ (A(2 n x) - A(2 n y)) for all x,y G M 2 . 

(i) Let us state the upper estimates first. We fix x,y G [0, l] 2 such that ||x — y\\ G (0, 1) and 
consider N ^ 1 such that 2~ N < \\x — y\\ ^ 2^( iV ~ 1 ). Using the fact that A is Lipschitz, with 
Lipschitz constant 1, and bounded by 1, we have \Z n (x,y)\ ^ min(||x — y\\, 1). Therefore 

oo N—l oo 

1/0*0 - f(y)\ < E y)\<T, ^ n(1 ~ a) i' x - yii) + E 2 ~ na < °W x - y\\ a 

n=0 n=0 n=N 

where C is a positive constant which only depends on a. Then, \f(x) — f(y)\ ^ Cr a for all 
x, y G M 2 such that < \\x — y\\ < \[2t < 1. This gives the upper bound in (|38|) . 

(ii) Now we state the lower estimates. We begin with finding a lower bound for the oscillation 
over an hexagonal cell. The key-point for estimating this oscillation is to calculate two 
particular increments for two pairs of well-chosen points belonging to the cell, namely the 
center and vertices. To keep in mind the number of the generation, we denote by C™, c G X n , 
the cells of the Voronoi tessellation of generation n. Moreover, the six vertices of a cell C" 
are denoted by Cj with i G { 1 , . . . , 6} . 
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Then, let N ^ 1, c € Xn and £ £ {1, . . . , 6}. We can write 



N-x 



/( c ) - /(ci) = X] ^ n ( c ' c ») + X] Zn ( c ' Q 



n=0 71=^ 

7JV 



As soon as a point is the center (resp. a vertex) of a cell then it is the center (resp. 
a vertex) of all the cells of higher generations. Hence A(2 n c) = 1 and A(2 n Cj) = for all 
n ^ N. Therefore 



oo 



n=N n=N 



z n (c,ci) = 2 

Let £, j £ {1, . . . , 6}. We obtain 

sup 1/(2/0 - /(y)| £ i((/(c) - /(<*)) + (/(c) - /(c,))) 



2 ( Z ^( c ' c i) + z »( c > c i) J+7" 



2 . 

n=0 n=0 

'A' 



-iVa 



For a fixed hexagon , we claim that there exists a pair of two diametrically opposite vertices 
Cj and Cj such that the sum of the two first sums above is positive. Indeed, for n < N, the 
center c is included in one or two cells C'J 1 of the tessellation of generation n and it can be 
only in three positions: at the center d , on one 'edge' (i.e. a segment between two consecutive 
vertices d i and c' i+1 or between the center d and its vertex c£), or on a 'face' (i.e. an open 
triangle of vertices d, d i and d i+1 ). Let us denote by C, E and F respectively these three 
positions. The sequence of positions when n goes from N — 1 to has to be (C, • • • ,C) or 
(C, • • • , C, E, • • ■ ,E)ot (C, C, • • • ,C,E, - ■ ■ ,E,F,--- , F) (see the first generations on Figure 
E] below). 




Figure 7. The only three possibilities for the sequence of the positions of the centers of 
the cells. On the left p goes from 1 to 0: c gives (C, C), c' gives (C, E) and c" gives (C, £?). 
On the right p goes from 2 to 0: c gives (C, C, C), c' gives (C, C, E), c" gives (C, -E, E) and 
c'" gives (C,E,F). 



In the first case, any pair of vertices will make the sum of the two sums above non-negative. 
In the second and third cases, we consider the uq when c is in position E. then, there are 
two diametrically opposite vertices associated with c on the same edge of the tessellation of 
order no and these two vertices satisfy the positivity of the sums above. 
Therefore, 

sup \F(y')-F(y)\ > C'2~ N <* 
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where C is a positive constant which only depends on a. 

Finally, let r G (0, 1), x G [0, l] 2 and N ^ 1 such that 2 _JV < r 2~( Ar - 1 ). The ball B T (x) 
contains a cell C^ +1 of generation TV + 1, thus 

osc T (x) > sup \F{y') - F(y)\ > c'2~ {N+ ^ a ^ C"r a , 
which gives the lower bound in (|38p . 

(ii) Combining (i) and Proposition 11.1 in [14] (see also [13]) we get A^(tat) ~ t^~ 3 as A" — >• oo. 
The result is then a consequence of (134j) . □ 



3.2. Takagi-Knopp series generated by a random perturbation of the dyadic mesh. 

In this subsection, we stray from the Voronoi partition of M. D . An alternative way of ran- 
domizing the underlying partition of a Takagi-Knopp series is the following: the sequence of 
dyadic meshes T> n = 2~ n lP , n G N, is kept but each mesh T> n is translated by a random 
uniform vector in 2~ n (0, 1) D . In each cube C U: k, k G N, of T> n , a random uniform 'nucleus' 
c n ^ is chosen independently. The associated random pyramidal function A n is defined so 
that it is equal to zero on T> n and to 1 on the set {c n ^,k G N}. The Takagi-Knopp type 
series F a ^, a, f3 > 0, is then given by ([2]). 

Compared to our Poisson- Voronoi construction, the main advantage of this model is that 
it preserves the cube structure so that it could be easier to deal with it in dimension two 
for applicational purpose in image analysis and pixel representation. The essential drawback 
is that the rigid structure of the mesh prevents us from obtaining an explicit Hausdorff 
dimension. Still, some results close to those proved in Section Q] can be deduced from similar 
methods. Indeed, if the random oscillation set B| jv, n,N G N, the increment Z n (x,y), 
x, y G R-°, the density z n of Z n (x, y) conditionally on {x G O n ^} and the Lipschitz constant 
L{x) are defined analogously, then the conclusions of Propositions II. ll and 11.31 are satisfied. 
Moreover, the point (ii) of Proposition 11.21 is replaced by the following estimate: 

„ D 
sup z n (t) — — - A — > -. 40 

te K nV P(x G O n , n ) ^ \ Xi - yi \ K ' 

As a consequence, it is possible to derive an analogue of Proposition 11.41 and of the upper- 
bound of the box-dimension in ([3]). This upper-bound is in particular the exact Hausdorff 
dimension when D = 1. Indeed, in the linear case, the estimate of sup t6R z n (t) in (|40p 
coincides with (ii) of Proposition 11.21 which implies that the lower-bound can be obtained 
along the same lines as in Section [2] We sum up our results in the next proposition, given 
without a detailed proof. 

Proposition 3.2. Let F a R be the function as above with < a ^ /3 ^ 1. Then, its graph 
r«,/3 satisfies almost surely the following estimates. 

(i) For every D > 2, dim B {T a ^) ^ D + 1 - |. 

(ii) When D = 1, dim B (r Qi/3 ) = dim H (r ai/3 ) = 2 - f . 

Unfortunately, the fact that the sum S^P , -, — - — r is a substitute for the inverse of the Eu- 

J ' L^i=\ \xi-yi\ 

clidean norm ||x — in the estimate (|40p makes the proof of the lower-bound of the 

Hausdorff dimension more intricate for D 2. This could be eventually considered as an 
extra-argument in favor of our Poisson- Voronoi construction. 
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